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Abstract.  Future  spaceborne  interferometric  arrays  must  meet  stringent 
optical  performance  and  tolerance  requirements  while  exhibiting  modu¬ 
larity  and  acceptable  manufacture  and  integration  cost  levels.  The  Mas¬ 
sachusetts  Institute  of  Teohnology  (MIT)  Adaptive  Reconnaissance 
Golay-3  Optical  Satellite  (ARGOS)  is  a  wide-angle  Fizeau  interferometer 
spacecraft  testbed  designed  to  address  these  research  challenges.  De¬ 
signing  a  space-based  stellar  interferometer,  which  requires  tight  toler¬ 
ances  on  pointing  and  alignment  for  its  apertures,  presents  unique  mul¬ 
tidisciplinary  challenges  in  the  areas  of  structural  dynamics,  controls, 
and  multiaperture  phasing  active  optics.  In  meeting  these  challenges, 
emphasis  is  placed  on  modularity  in  spacecraft  subsystems  and  optics 
as  a  means  of  enabling  expandability  and  upgradeability.  A  rigorous 
theory  of  beam-combining  errors  for  sparse  optical  arrays  is  derived  and 
flown  down  to  the  design  of  various  subsystems.  A  detailed  elaboration 
on  the  optics  system  and  control  system  is  presented  based  on  the  per¬ 
formance  requirements  and  beam-combining  error  tolerances.  The 
space  environment  is  simulated  by  floating  ARGOS  on  a  frictionless  air¬ 
bearing  that  enables  it  to  track  both  fast  and  slow  moving  targets.  ©  2004 
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1  Introduction 

1.1  Sparse  Aperture  Interferometric  Arrays 

The  quest  for  finer  angular  resolution  in  astronomy  inevi¬ 
tably  leads  to  larger  aperture  since  the  resolution  is  propor¬ 
tional  to  the  wavelength  over  the  diameter  of  a  circular 
aperture.  Unfortunately,  the  primary  mirror  diameter  for 
space  telescopes  is  limited  by  volume  and  mass  constraints 
of  current  launch  vehicles  (ca.  4  to  5  m)  as  well  as  the 
scaling  laws  of  manufacturing  costs.*  Since  the  cost  of 
monolithic  optics  increases  faster  than  their  diameter 
squared,  and  mirrors  such  as  the  Hubble  Space  Telescope’s 
are  already  at  the  edge  of  what  is  financially  feasible,  ef¬ 
forts  are  ongoing  to  break  this  trend  by  employing  break¬ 
through  technologies  such  as  deployed  segmented  mirror 
telescopes  and  sparse  aperture  optics  using  interferometry. 
Whereas  long  baseline  stellar  Michelson^  interferometers 
feed  light  from  independent  collectors  to  a  beam  combiner 
to  obtain  interference  fringes  over  a  period  of  time  (pupil 
plane  beam  combination),  Fizeau^  interferometers  produce 
a  direct  image  with  full  instant  u-v  coverage  (image-plane 
beam  combination).  Hence,  Fizeau  interferometers  are  suit¬ 
able  for  optical  imaging  of  extended  objects  and  rapidly 
changing  targets.  In  contrast  to  the  long  baselines  of  Mich- 
elson  interferometers,  Fizeau  interferometry  systems  tend 
to  have  compact  telescope  arrays.  An  optimal  imaging  con¬ 


figuration  designed  for  sparse  arrays  was  first  proposed  by 
Golay.^  Sparse  arrays  are  promising  for  applications  that  do 
not  require  extremely  high  sensitivity  (bright  source 
present)  and  allow  for  a  rather  limited  field  of  view"*’^ 
(FOV).  Diffraction-limited  performance  has  been  demon¬ 
strated  for  sparse  or  dilute -aperture  telescopes  with  active 
phasing  control.®  A  notable  project  in  the  area  of  phased 
telescope  arrays  is  the  Multipurpose  Multiple  Telescope 
Testbed’  (MMTT)  built  by  the  Air  Force  Research  Labora¬ 
tory  (AFRL).  The  MMTT  consists  of  four  20-cm-aperture 
telescopes  phased  together  with  a  15-arcmin  FOV.  The 
MMTT  employs  a  complex  laser  interferometer  metrology 
to  sense  wavefront  error  (WFE).  The  Multi-Aperture  Imag¬ 
ing  Array*  built  by  Lockheed  Martin  demonstrated  phase 
diversity  computation  techniques  for  WFE  sensing.  This 
sparse  array  consists  of  afocal  telescopes  arranged  in  a  Y 
formation  that  are  combined  to  a  common  focus  in  a  Eizeau 
interferometer  configuration.  It  demonstrated  the  first  re¬ 
sults  of  a  broadband  multiple-telescope  imaging  array 
phased  over  a  significant  EOV  using  an  extended  image 
projector  in  the  lab.  Research  into  WEE  sensing  and  control 
has  been  extensively  conducted  for  the  Next  Generation 
Space  Telescope  (renamed  the  James  Webb  Space 
Telescope).® 

The  Massachusetts  Institute  of  Technology  (MIT)  Space 
Systems  Laboratory  (SSL)  successfully  completed  the  Mid- 
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Fig.  1  Golay-3  ARGOS  system  with  the  three  attitude  control  sys¬ 
tem  (ACSs)  shown  at  the  bottom. 
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deck  Active  Control  Experiment**'  (MACE),  a  Space 
Shuttle  flight  experiment  that  flew  on  STS-67  in  March 
1995.  The  primary  goal  of  the  MACE  project  was  to  inves¬ 
tigate  approaches  for  achieving  high-precision  pointing  and 
vibration  control  of  future  spacecraft.  The  MIT-SSL  also 
developed  The  Origins  Testbed**  (OT),  a  laboratory  test 
article  to  capture  the  dynamics  and  control  problems  antici¬ 
pated  for  future  space  telescopes.  Built  on  the  heritage  of 
the  MACE  and  the  OT,  a  ground-based  satellite  testbed 
with  a  more  sophisticated  active  optical  payload  system  is 
presented  in  this  paper.  Most  of  the  previous  sparse- 
aperture  systems  do  not  deal  with  real-world  problems  such 
as  the  vibrational  coupling  between  a  spacecraft  structure 
and  the  wavefront  errors  propagating  through  the  whole 
system.  The  presented  testbed  explores  such  issues  while 
implemented  on  a  nonrotationally  constrained  platform. 

1 .2  Overview  of  ARGOS 

To  better  understand  the  technological  difficulties  involved 
in  designing  and  building  a  sparse  aperture  array,  the  chal¬ 
lenge  of  building  a  Golay-3  interferometer  system  was  un¬ 
dertaken.  The  MIT  Adaptive  Reconnaissance  Golay-3  Op¬ 
tical  Satellite*^  (ARGOS)  exploits  wide-angle  Eizeau 
interferometer  technology  with  an  emphasis  on  modularity 
in  the  optics  and  spacecraft  subsystems.  The  objective  of 
the  ARGOS  project  is  to  demonstrate  the  feasibility  of  a 
modular  architecture  for  space-based  optical  systems. 

To  demonstrate  a  complete  spacecraft  in  a  1  -,g  environ¬ 
ment,  the  ARGOS  system  is  mounted  on  a  frictionless  air¬ 
bearing,  as  shown  in  Eig.  1,  and  has  the  ability  to  track  fast 
orbiting  satellites  like  the  International  Space  Station  (ISS) 
as  well  as  slow  moving  objects  such  as  point  stars.  The 
modular  architecture  emphasizes  the  use  of  replicated  com¬ 
ponents  and  standard  interfaces.  The  system  consists  of 
three  identical  apertures  arranged  in  a  Golay-3  distribution. 
The  light  from  these  telescopes  is  combined  in  a  center 
module  and  transmitted  to  a  charge-coupled  device  (CCD). 


Table  1  Key  functional  requirements. 


Key  Requirements 

Angular  resolution 

0.35  arcsec  at  operating  wavelength 

Operating  wavelength 

400-700  nm  (visible) 

FOV 

3  arcminx3  aremin 

Field  of  regard  (FOR) 

120  deg 

Signal-to-noise  ratio  (SNR) 

100 

ACS  pointing  accuracy 

±  1  aremin 

Image  acquisition  time 

20  images/h  (max) 

Autonomous  operation  time 

up  to  1  continuous  h 

Wavefront  sensing  techniques  are  explored  to  mitigate  ini¬ 
tial  misalignment  and  to  feed  back  real-time  aberrations 
into  the  optical  control  loop.  The  goal  is  to  obtain  an  image 
of  similar  quality  as  the  image  received  from  a  monolithic 
telescope  using  a  single  comparable  aperture.  The  primary 
functional  requirements  of  ARGOS  are  prescribed  in  Table 
1. 

The  structure  of  this  paper  follows  the  step-by-step  de¬ 
sign  procedures  of  ARGOS.  A  theory  of  interferometric 
beam  combining  errors  is  developed  in  the  process  and 
connected  to  the  actual  design  and  implementation  of  AR¬ 
GOS.  The  paper  also  details  the  design  of  the  optics  based 
on  the  performance  requirements  and  beam-combining  tol¬ 
erances.  In  addition,  it  includes  the  control  avionics  and 
strategies  for  a  large-angle  slewing  and  active  interferomet¬ 
ric  beam  combining. 

2  Physics  of  Sparse-Aperture  Interferometric 
Arrays 

2.1  Determination  of  Array  Configuration 

Traditional  image  quality  criteria  such  as  resolution  and 
encircled  energy  (EE)  are  inadequate  for  many  sparse- 
aperture  or  interferometric  array  applications.*^  When  look¬ 
ing  at  extended  objects  such  as  the  Moon  or  faint  distant 
nebulae,  evaluation  of  an  optical  system  involves  more  than 
simply  looking  at  a  point  source  response  [point  spread 
function  (PSE)].  The  modulation  transfer  function  (MTE)  is 
a  better  metric  to  evaluate  the  contrast  (modulation)  trans¬ 
fer  characteristic  of  an  extended  object.  Eigure  2  shows  the 
crucial  relationship  between  the  two  most  common  image 
quality  criteria — PSE  and  MTE.  The  PSE  is  the  squared 
modulus  of  the  Eourier  transform  of  the  complex  pupil 
function.  The  optical  transform  function  (GTE)  is  a  Eourier 
transform  of  the  PSE,  and  the  MTE  is  an  absolute  value 
(magnitude)  of  the  GTE.  In  Eig.  3,  the  PSE  and  MTE  plots 
of  a  GoIay-3-like  array  with  D  =  0.2I  m  (ARGGS  subaper¬ 
ture  diameter)  and  the  array  radius,  L  =  0. 12  m  (solid),  0.19 
m  (dash),  0.3  m  (dot)  are  shown.  A  perfect  monolithic  ap¬ 
erture,  free  of  optical  aberrations,  has  a  linearly  decreasing 
MTE  contrast  characteristic  (see  the  solid  MTE  line  in  Eig. 
3).  In  case  of  a  sparse  array,  the  MTE  suffers  a  contrast  loss 
in  the  mid  spatial  frequency  range,  as  shown  in  Eig.  3.  We 
can  observe  that  the  MTE  plot  with  L  =  0.3  m  exhibits  two 
zero  MTE  values  rather  than  one.  The  first  zero  denoted  by 
F ^  is  the  practical  spatial  cutoff  frequency,  and  defines  the 
“practical  resolution  limit.”  The  is  the  cutoff  frequency. 
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Fig.  2  The  relationship  between  the  PSF  and  OTF. 


whose  inverse  indicates  an  angular  resolution  under  the 
normal  condition  that  there  is  no  F ^  (another  zero  point) 
before  .  Thus,  the  larger  F^  or  F ^  is,  the  better  is  the 
angular  resolution  a  sparse  array  will  achieve. 

In  contrast  with  monolithic  apertures,  the  method  of  full 
width  at  half  maximum  (FWHM)  is  not  sufficient  to  deter¬ 
mine  the  angular  resolution.  Assuming  that  angular  resolu¬ 
tion  is  fully  determined  by  the  array  size,  the  PSF  can  re¬ 
veal  the  highest  achievable  angular  resolution.  This 
assumption  holds  especially  for  very  large  baseline  Mich- 
elson  interferometers.  Figure  3  shows  that  the  width  of  the 
mainlobe  of  the  PSF  plot  is  getting  smaller,  indicating  im¬ 
proving  angular  resolution  as  we  increase  the  array  size  L. 
The  Fizeau  interferometer,  however,  requires  an  instanta¬ 
neous  full  u-v  coverage,  which  limits  the  practical  resolu¬ 
tion.  In  addition,  as  discussed,  considerable  contrast  loss  in 
the  mid-spatial-frequency  range  should  be  avoided.  The  an¬ 
gular  resolution  from  the  MTF  plot  of  a  Golay-3  type 
sparse  imaging  array  is  approximated  as*^ 


\  \ 

=  1  22 - =  2  44 - , - ^ 

'■  D,s  3L  +  (4D2-3L2)1/2' 


(1) 


As  L  is  increased,  the  array  becomes  sparser,  which,  as  a 
result,  boosts  the  heights  of  sidelobes  in  the  PSF.  The  MTF 
plot  at  L  =  Q3  m  (dotted  line  in  Fig.  3)  has  two  zeros  while 
others  have  only  one.  In  that  case,  the  resolution  is  limited 
not  by  F^  but  by  F,. ,  whose  inverse  is  the  angular  resolu¬ 
tion  of  a  single  aperture  (1.22  k/D).  Thus,  in  that  case,  the 
sparse  array  has  no  advantage  over  a  single  subtelescope 
except  for  increased  light-gathering  power.  The  array  con¬ 
figuration  of  ARGOS  is  selected  as  L  =  0.191 85  m  for  D 
=  0.21  m  subtelescopes.  This  results  in  a  better  theoretical 
angular  resolution  of  0.35  arcsec  for  the  550-nm  wave¬ 
length  rather  than  0.55  arcsec  of  a  single  aperture,  compro¬ 
mising  with  a  reasonable  MTF  characteristic. 


2.2  Beam  Combining  Errors 

There  are  three  major  wavefront  errors  that  must  be  con¬ 
trolled  at  the  beam  combiner  to  achieve  phased  beam  com¬ 
bining.  Those  errors  are  optical  path  difference  (OPD),  tip/ 
tilt  error,  and  pupil-mapping  error. 


PSF  Plot 


Fig.  3  PSF  and  MTF  plots  when  D=0.21  m  and  L=0.12  m  (solid),  0.19  m  (dash),  and  0.3  m  (dot). 
The  corresponding  array  configurations  are  shown  to  the  right  of  PSF-MTF  plot.  The  black-gray  figure 
is  an  MTF  plot.  Gray  circles  indicate  the  practical  cutoff  frequency. 
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(a)  (b)  (c) 

Fig.  4  PSF  plot  of  Golay-3  array  with  (a)  0  OPD,  (b)  0.5\  OPD,  and 
(c)  1.0\  OPD. 


2.2.1  OPD  (piston)  error 

We  can  plot  the  effects  of  OPD  errors  using  the  interferom¬ 
etry  equation  given  by  Mennesson  and  Mariotti*"^ 


1 - 1 

o 

o 

+ 

J lirrD  sin  r/X) 

X 

n 

ttD  sin  r/X 

X 


2  exp[y2'7r(Li.r/\)cos(<5i.— 
k=l 


2 


(2) 


where  r  is  the  off-axis  angular  direction,  0  is  the  azimuth 
angle,  D  is  a  subaperture  diameter,  {Li^,Si()  are  the  polar 
coordinates  of  the  array  configuration,  and  is  the  phase 
shift. 

Numerical  simulations  using  Eq.  (2)  were  performed  to 
investigate  the  effects  of  OPD  errors,  as  shown  in  Fig.  4.  As 
the  piston  error  increases,  two  major  deviations  develop 
over  the  envelope  of  the  PSF.  First,  the  main  envelope 
shifts  in  the  direction  of  the  piston  error.  The  resultant  di¬ 
rection  of  the  envelope  shift  is  the  vector  sum  of  piston 
error  directions  weighted  by  the  amount  of  error.  Second, 
the  peak  intensity  is  reduced  compared  to  the  nominal  PSF 
without  any  piston  errors  resulting  in  a  reduced  Strehl  ratio 
(SR).  The  size  of  the  mainlobe  also  expands  showing  a 
degraded  angular  resolution  [see  Fig.  4(c)].  When  the  pis¬ 
ton  error  is  O.IX.  between  apertures,  the  peak  intensity  is 
98%  of  the  normal  intensity.  The  beam  combining  piston 
error  tolerance  is  chosen  to  be  O.IX  =  55  nm  for  most  of 
interferometric  beam  combining  applications. 

2.2.2  Tip/tilt  error 

The  approach  employed  here  to  analyze  tilt  errors  is  to 
further  segment  each  aperture  into  smaller  elements  for 


analysis  (finite  element  method).  The  PSF  is  calculated  by 
summing  up  the  interference  from  each  element.  The  phase 
difference  at  the  central  point  due  to  tilt  errors  are  added  to 
the  interference  term  of  Eq.  (2): 


T7C/[1  -l-cos(r)] 

J lirrd  sin  r/X) 
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Trd  sin  r/X 
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2  exp[;27r(L,tr/X)cos(5^-0)] 
k=l 
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Xexpljdk  sin(tilt„)27r/X] 


(3) 


where  D  in  Eq.  (2)  is  replaced  with  the  smaller  element 
diameter  d,  d/.  is  the  distance  from  the  aperture  center  to 
the  k’th  element  center,  tilt,,  denotes  the  tilt  in  radians  for 
the  «’th  aperture,  k  ranges  from  1  to  m,  and  m  equals  the 
number  of  apertures  («)  multiplied  by  the  number  of 
smaller  elements  of  each  aperture.  By  reading  the  SR  val¬ 
ues  from  the  PSF  plots  under  the  influence  of  tilt  errors,  the 
maximum  allowable  tilt  error  at  the  beam-combining  stage 
is  determined.  This  numerical  method  predicts  that  the  tilt 
tolerance  between  each  beam  entering  the  beam  combiner 
should  be  less  than  20  /rdeg  (0.072  arcsec=0.35  /rrad). 


2.2.3  Pupil  mapping  error 

If  coherent  imaging  is  to  be  achieved  over  any  significant 
FOV,  the  pupil  mapping  process  must  be  performed  such 
that  the  exit  pupil  is  an  exact  (scaled)  replica  of  the  en¬ 
trance  pupik^  This  constraint  is  commonly  called  the 
“golden  rule”  of  beam  combining.  From  the  geometry  of 
Fig.  5,  Faucherre  et  al.^  derived  the  net  OPD  error  due  to 
the  incorrect  pupil  mapping. 


OPD„g(=  \B  sina  —  b  sin  /3\ 


B 

B  sin  a - sin(m„a') 


(4) 


where  is  the  aperture  magnification  factor  (Did),  is 
the  baseline  magnification  factor  {Bib),  a  is  half  FOV  of 
the  target,  and  fd  is  magnified  field  angle  through  the  sub¬ 
telescope  array. 
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Fig.  5  The  golden  rule  of  beam  combining,  pupil  mapping. 
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According  to  the  discussion  in  Sec.  2.2.1,  1/lOX.  phasing 
is  desired  to  meet  the  SR  and  angular  resolution  require¬ 
ments.  Thus,  equating  Eq.  (4)  with  1/10\  and  small  angle 
approximation  sin  a=a  results  in 


Ba 
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1 - ^ 
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\ 

To' 


(5) 


A  calculation  of  shear  error  tolerance  is  performed  with  an 
assumption  that  there  is  no  magnification  error  (subtele¬ 
scopes  are  strictly  afocal  with  exactly  the  same  magnifica¬ 
tion).  Then,  the  correct  magnification  is  m^  =  D/d  =  B/(b 
+  Ab).  The  incorrect  mapping  is  still  Bib.  Inserting 
and  rris  into  Eq.  (5)  results  in 
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Fig.  6  FOV  versus  pupil  mapping  tolerances  for  ARGOS  with 
achievable  magnification  tolerance. 


Solving  for  Ab,  gives 


Ab  = 


X 

lOm^a 


(7) 


Equation  (7)  means  that  the  allowable  shear  error  is  in¬ 
versely  proportional  to  both  aperture  magnification  and 
off-axis  target  angle  a  (EO'V/2).  However,  the  assumption 
of  no  magnification  error  used  for  Eq.  (7)  is  inadequate 
since  perfectly  replicating  any  subaperture  is  impossible  in 
terms  of  manufacturability.  There  usually  exists  some  per¬ 
formance  window  specified  by  optics  manufacturers.  The 
objective  of  this  calculation  is  to  find  out  how  precise  the 
magnification  tolerance  level  must  be  in  order  to  achieve 
coherent  beam  combining.  If  the  magnification  tolerance  is 
too  stringent  to  be  manufactured,  methods  of  relaxing  this 
tolerance  will  be  explored. 

Eor  both  incorrect  m^{  =  Dld)  and  incorrect  m, 
{  =  B/b),  the  acceptable  magnification  ratios  are  defined  as 
follows: 
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Inserting  these  results  into  Eq.  (5),  gives 


ab 


B 

B  B 

- 

=  ab 

- 7 - h  Am„ 

b 

b  b  +  Ab 

AbB 


-  +  Am„b 


b  -F  Ab 

We  can  remove  b  by  using  Eq.  (8): 
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Setting  q;  =  EOV/2,  and  B  =  V2L  and  for  a 

Golay-3  array  leads  to 
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(11) 


Erom  Eq.  (11),  we  can  find  the  EO'V  as  a  function  of  mag¬ 
nification  (m)  and  pupil  mapping  tolerances  (A/  and  Am): 


EOV= - .  (12) 

5V2\mAl+  {AmL/m)  —AmAl\ 

We  can  also  represent  the  shear  error  tolerance  A1  as  a 
function  of  EOV,  magnification  (m)  and  magnification  tol¬ 
erance  (Am)  for  a  positive  Am: 


X 

Al= - 

5v3EOV(m-Am) 


AmL 

m^  —  niAm  ' 


(13) 


2.2.4  Pupil  mapping  error  for  ARGOS 

These  results  are  graphically  plotted  in  Eigs.  6  and  7.  At  the 
preliminary  design  phase  of  ARGOS,  12  /jm  was  suggested 
for  the  shear  error  tolerance  from  Eq.  (7)  with  an  assump¬ 
tion  of  no  magnification  error  (EOV  =  4  aremin, 
magnification  =  10).  However,  12-yU,m  shear  error  cannot 
produce  the  EOV  requirement  of  ARGOS  (3  aremin)  for 
any  range  of  magnification  error  in  Eig.  6.  We  can  tighten 
shear  error  tolerance  to  meet  the  EOV  requirement  or  we 
can  relax  the  EOV  requirement  by  shrinking  the  region  of 
interest  within  the  whole  EOV.  Therefore,  the  pupil  map¬ 
ping  process  is  the  primary  limiting  factor  deciding  the  rea¬ 
sonable  EOV  of  a  sparse  aperture  interferometric  array.  AR¬ 
GOS ’s  subtelescope  collimators  are  designed  with  a 
tolerance  of  ±0.0095  with  a  nominal  magnification  factor 
of  10.  However,  the  ARGOS  subtelescopes  have  a  focusing 
knob  that  can  control  the  distance  between  the  primary  mir¬ 
ror  and  the  secondary  mirror,  thereby  controlling  the  size  of 
the  beam  more  precisely.  It  is  usually  considered  that  a 
1/1000  magnification  tolerance  requirement  is  too  expen¬ 
sive  to  manufacture.  Eigure  7  describes  how  the  magnifica¬ 
tion  of  a  subaperture  affects  the  shearing  tolerance  of  a 
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Magnification  tolerance  l-Oe-S 


^  ,  ,  .  Magnification  of  Subaperture 

Shear  Error  Tolerance  [microns] 

Fig.  7  Magnification  versus  shear  tolerance  with  magnification 
error=  0.001. 


sparse  aperture  imaging  system.  This  plot  implies  that  sub¬ 
aperture  magnification  can  be  tuned  to  maximize  allowable 
shear  error  A/  (lateral  pupil  mapping  error)  thereby  reduc¬ 
ing  control  complexity.  For  a  magnification  of  10,  we  can 
increase  the  shear  tolerance  value  sufficiently  high  at  the 
expense  of  the  reduced  FOV. 

2.3  Optics  Design  and  Impiementation 

2.3.1  Subaperture 

The  overall  cost  of  the  ARGOS  optics  system  could  be 
significantly  reduced  by  selecting  one  of  the  highest  preci¬ 
sion  optics  commercially  off-the-shelf  (COTS)  telescope. 
However,  it  was  necessary  to  customize  a  collimating  lens 
to  convert  the  Dall-Kirkham-type  focal  telescope  to  an  afo- 
cal  telescope  with  a  magnification  ratio  of  10.  The  collimat¬ 
ing  lens  is  placed  into  the  baffle  of  the  telescope  to  make 
the  system  compact.  A  trade-off  analysis  was  performed  on 
several  different  cemented  doublets  that  were  optimized  ex¬ 
tensively  by  ZEMAX.  One  drawback  of  the  cemented  dou¬ 
blet  is  that  it  has  bonded  glasses,  therefore  if  there  is  a 
change  of  temperature,  the  doublet  may  fail.  Although  a 
doublet  with  CaF2  performs  best  in  reducing  chromatic  ab¬ 
errations,  the  high  coefficients  of  thermal  expansion  (CTEs) 
of  CaE2  (18.3)  forced  to  find  another  glass  combination  for 
efficient  achromatic  doublet  design.  Smith^^  suggests  EK51 
(as  a  crown  element)  with  a  KzES  or  LaK  glass  (as  a  flint). 
Although  the  maximum  focal  shift  range  can  be  minimized 
to  247  yLtm  with  EK51-KzESll,  it  was  not  the  best  choice 
due  to  the  residual  aberrations  [root  mean  square  (rms) 
wavefront  errors  predicted  by  ZEMAX].  The  final  EK51- 
BaK2  design  achieves  271.6  pm  chromatic  focal  shift 
range. 

2.3.2  Design  of  OPD  and  tiit/tip  controiier 

When  strictly  looking  at  two  designs  of  optical  delay  lines 
(ODE)  in  Eig.  8,  a  perpendicular  design  seems  to  offer  the 
most  benefits.  Eor  the  perpendicular  design,  a  multiaxis  fast 


Fig.  8  Perpendicular  ODL  design  with  FSM  (left)  and  parallel  ODL 
design  coupled  with  FSM  (right). 

steering  mirror  (ESM)  could  be  used  to  generate  the  re¬ 
quired  tip/tilt  actuations  as  well  as  fine  optical  path  differ¬ 
ence  control.  ESMs  have  very  fine  resolutions  so  depending 
on  them  for  fine  OPD  control  will  enable  us  to  have  a  more 
reasonable  coarse  control  for  the  ODL.  When  the  parallel 
ODL  design  is  coupled  with  a  ESM  (Eig.  8),  the  resulting 
design  is  very  simple  and  more  cost-effective  than  the  per¬ 
pendicular  design.  This  design  cannot  perform  coarse  OPD 
control,  but  this  design  is  more  compact  resulting  in  easier 
integration  with  the  structural  design.  In  addition,  there  is 
no  need  for  a  translational  stage  or  rooftop  mirror,  which 
reduces  the  cost  and  control  complexity,  and  there  is  also  a 
greater  total  reflectance  since  there  are  fewer  mirrored  sur¬ 
faces.  Eewer  mirrored  surfaces  also  lead  to  fewer  structural 
misalignment  errors.  In  both  the  perpendicular  ODL  design 
and  the  parallel  ODL  design,  fine  OPD  control  is  coupled 
with  shear  control  so  that  AEine  OPD  =  AShear.  The  effect 
that  a  change  in  fine  OPD  would  have  on  shear,  is  not  very 
significant  and  could  be  ignored  for  adjustments  <10  pm. 
This  is  because  we  have  a  much  tighter  tolerance  on  piston 
error  (50  nm)  than  the  \2-pm  shear  error.  A  fine  resolution 
multiaxis  ESM,  capable  of  controlling  the  tilt/tip  as  well  as 
piston  motion,  was  selected  as  the  ODL-ESM  actuator  of 
ARGOS.  The  ESM  has  to  be  able  to  compensate  for  any 
errors  in  its  mounting.  Therefore  a  high-precision  mount 
with  a  range  up  to  7  deg  was  selected  for  the  ESM  actuators 
(see  Eig.  10  and  Table  2  in  Sec.  2.3.3). 

2.3.3  Pyramidai  mirror  and  beam  combiner 

The  pyramidal  mirror  turns  all  three  beams  45  deg  into  the 
beam  combiner.  A  custom  pyramidal  mirror  is  chosen  due 
to  the  cost  of  making  one  out  of  regular  mirrors.  The  main 
reason  for  the  high  cost  is  thin  mirrors  that  cost  $10003- , 
and  then  mount  them  to  an  accuracy  of  ±0.001  deg 
( ±  3.6  arcsec).  The  pyramid  cannot  be  made  from  the  regu¬ 
lar  thickness  mirrors  since  they  constrain  the  beam  diam¬ 
eter.  We  customized  a  pyramidal  mirror  with  a  surface  ac¬ 
curacy  of  \/10  peak  to  valley,  ±3  arcsec  angle  error,  and  a 
50-mm  clear  aperture.  The  substrate  material  is  BK7  with  a 
coating  of  AlSiO  (aluminum  with  silicon  monoxide).  The 
reflectance  will  be  approximately  90%  in  the  visible  range. 
The  two  PSEs  shown  in  Eig.  9  demonstrate  how  the  ESM 
can  compensate  for  the  ±  3  arcsec  errors  in  the  pyramid.  To 
compensate  for  pyramid  errors  the  ESM  has  to  align  itself 
so  that  the  two  reflecting  surfaces  are  parallel.  The  pyrami¬ 
dal  mount  is  composed  of  two  stages.  The  first  stage  pro¬ 
vides  all  of  the  angular  adjustments  and  the  second  handles 
X  and  Y  translation  in  the  entrance  pupil  of  the  beam  com¬ 
biner.  The  X-Y  translation  stage  is  small  enough  to  fit  be¬ 
hind  the  tip/tilt  rotation  stage  and  has  the  load  capacity  to 
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Fig.  9  PSF  simulation  of  ARGOS  with  3-arcsec  tilt  error  of  the  py¬ 
ramidal  mirror  (SR=  0.444,  left).  The  PSF  with  FSM  correction 
(SR=  0.960,  right). 


hold  both  the  second  stage  and  the  mirror.  Any  additional  Z 
directional  error  can  be  offset  by  the  FSM  mounts.  The 
specifications  of  the  ARGOS  optical  actuators  and  mounts 
are  tabulated  in  Table  2  and  their  pictures  are  shown  in  Fig. 
10. 

Two  options  available  for  the  image-plane  beam  com¬ 
biner  are  either  reflecting  or  refracting  optics.  The  reflecting 
beam  combiner  is  compact  when  compared  to  a  refractor. 
Unfortunately,  the  secondary  mirror  of  a  Cassegrain  tele¬ 
scope  would  partially  block  the  three  incoming  beams  in 
any  possible  configurations  to  obey  the  golden  rule  dis¬ 
cussed  in  Sec.  2.2.3.  A  single  parabolic  mirror  was  consid¬ 
ered,  however,  there  was  not  enough  space  between  the 
pyramidal  mirror.  Had  we  used  reflecting  optics,  there 
would  be  nothing  available  as  COTS,  therefore  it  would 
have  to  be  custom  manufactured  increasing  cost  signifi¬ 
cantly.  A  reflector  would  also  complicate  the  relay  optics 
significantly  since  we  would  not  be  able  to  use  the  pyrami¬ 
dal  mirror.  In  contrast,  a  refracting  combining  telescope  has 
many  advantages.  It  allows  for  very  simple  relay  optics.  It 
is  available  as  COTS  with  high-quality  optics  such  as  Ta- 
kahashi  FSQ-106N  (see  Fig.  10).  This  telescope  has  signifi¬ 
cantly  less  chromatic  aberration  than  other  COTS  tele¬ 
scopes.  This  is  due  to  its  four-element  design,  two  of  which 
are  fluorite.  It  has  a  diameter  of  106  mm,  and  a  530-mm 
focal  length,  resulting  in  a  total  system  focal  length  of  5300 
mm. 

2.3.4  CCD  system  design 

The  optics  performance  requirements  driving  the  CCD 
camera  selection  are  SNR  and  FOV.  It  was  pointed  out  that 
the  pupil  mapping  process  is  the  primary  limiting  factor 
deciding  the  reasonable  FOV  of  a  sparse  aperture  interfero¬ 
metric  array  (see  Sec.  2.2).  This  means  that  the  pupil  map¬ 
ping  tolerances  become  expensive  as  we  increase  the  FOV. 


Table  2  Optical  actuators  and  mounts  specifications. 


Model 

Angular  Range 

Angular 

Resolution 

Linear 

Range 

Linear 

Resolution 

FSM 

±600  yitrad 

±0.05  /irad 

12  /im 

0.2  nm 

FSM  mount 

±7  deg 

±0.0008  deg 
(±14  yiirad) 

1  cm 

1  fjm 

Pyramid  mirror 
mount 

±4  deg 

±2  arcsec 
(±9.6  Atrad) 

13  mm 

3  ijm 

Fig.  10  Pictures  of  optical  components. 


In  addition,  the  off-axis  imaging  aberrations  such  as  distor¬ 
tion  and  held  curvature  become  dominant  as  we  increase 
the  FOV.  Therefore,  the  CCD  system  of  a  sparse  aperture 
imaging  system  should  be  carefully  designed  to  meet  the 
optics  performance  requirements  while  maintaining  the  rea¬ 
sonable  FOV  large  enough  to  capture  the  targets  of  interest. 
The  following  relationship  is  used  to  And  an  appropriate 
CCD  pixel  size  and  the  number  of  CCD  pixels: 

d  \ 

¥OV=-n=2M  —  Qn,  (14) 

J  ^eff 

where  n  is  the  number  of  pixels  along  the  x  axis  or  the  y 
axis  of  a  CCD  matrix,  d  is  the  pixel  size,  /  is  the  system 
effective  focal  length,  and  Q  is  a  quality  factor.  The  Ny- 
quist  sampling  criterion  tells  us  that  two  pixels  are  gener¬ 
ally  required  to  properly  record  a  star’s  image  (2  =  0.5). 

A  sparse  aperture  can  cause  a  substantial  depression  in 
the  MTF  of  the  system,  resulting  in  a  low-contrast  and 
unsatisfactory  SNR  image.  This  can  be  compensated  by 
image  reconstruction,  such  as  by  using  a  Wiener-Helstron 
Alter,  but  at  the  expense  of  greater  noise  sensitivity  com¬ 
pared  to  a  monolithic  aperture.*®  Images  acquired  by  mod¬ 
ern  CCD  cameras  may  be  contaminated  by  a  variety  of 
noise  sources.  Three  primary  noise  sources  are  discussed  in 
this  section:  photon  noise,  read-out  noise,  dark  current. 

When  the  physical  signal  detected  by  a  CCD  is  based  on 
light,  the  quantum  nature  of  light  plays  a  significant  role. 
Photon  noise  is  derived  by  calculating  the  number  of  pho¬ 
tons  received  through  an  imaging  system,  and  is  therefore 
not  independent  of  light  source.  The  spectral  energy  distri¬ 
bution  of  a  blackbody  is  given  by  Planck’s  law:*^ 

1 

exp(c/7/k7’\)-l  ’ 

where  E{X)  is  the  energy  per  unit  wavelength  (also  called 
the  spectral  irradiance  measured  in  Wm^^  /rm^*),  X.  is  the 
wavelength,  h  is  Planck’s  constant  (6.6260755 
X  10^^^  W  s^),  T  is  the  absolute  temperature  of  a  target,  c 
is  the  speed  of  light,  and  k  is  Boltzmann’s  constant 
(1.380658X10^^^  Ws/K). 
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Then  the  input  power  at  the  entrance  to  a  sparse  aperture 
array  is  calculated  by 


P 


where  the  integration  term  integrates  E{\)  from  Eq.  (15) 
over  the  operating  wavelengths,  N  is  the  number  of  subap¬ 
ertures,  D  is  the  diameter  of  a  subaperture,  and  is  the 
effective  diameter  of  a  sparse  aperture  array. 

Using  Eq.  (16),  the  number  of  available  photons  can  be 
determined  as  follows: 

=  (17) 


iDYi\.22\Y 


(16) 


where  T,  is  the  integration  time,  Tq  is  the  optical  transmis¬ 
sion  factor  (typically  0.75),  c  is  the  speed  of  light,  and  h  is 
the  Planck  constant  dehned  in  Eq.  (15). 

The  number  of  signal  electrons  is  computed  by  multi¬ 
plying  with  the  CCD  quantum  efficiency: 

lVsignal=-Pin|  ^O^QEjT,. 

r>^2  I  D^\i\.22\Yl  \ 

^J^^E(X)dx(iv^J(— )  (.o^QE)r,, 

(18) 

where  QE  is  the  quantum  efficiency  of  a  CCD.  The  term 
N(D  /D^g)  is  commonly  called  the  hlling  factor  of  an  in¬ 
terferometric  array. 

The  probability  distribution  for  the  average  number  of 
photons  received  by  a  CCD  is  known  to  be  Poisson.  The 
standard  deviation  of  the  Poisson  distribution  is  the  square 
root  of  the  average  value.  Accordingly,  the  total  number  of 
photon  noise  electrons  is 


A^pn— 


signal 


£'(\)dX  N  — 


D 


eff, 


1.22\ 


-O^QEIT, 


1/2 


(19) 


By  cooling  the  CCD  chip  it  is  possible  to  efficiently  reduce 
the  number  of  thermal  electrons  that  produce  thermal  noise 
or  dark  current.  As  the  integration  time  increases,  the  num¬ 
ber  of  thermal  electrons  increases: 


1Vdc“^dcE,, 


(20) 


where  is  the  dark  current  rate  and  T,  is  the  integration 
time. 

Readout  noise  (on-chip  electronic  noise)  originates  in 
the  process  of  reading  the  signal  from  the  sensor  through 
the  held  effect  transistor  (PET)  of  a  CCD  chip.  This  noise 
can  become  a  signiheant  component  in  the  overall  SNR  at 
very  low  signal  levels. 

Pinally,  we  can  calculate  the  SNR  of  a  sparse  aperture 
array  based  on  the  noise  characteristics  of  a  CCD  system, 
as  already  discussed.  In  general,  photon  noise,  dark  current 
noise,  and  readout  noise  are  considered  as  dominant  noise 
contributors: 


SNR= 


N. 


signal 


(Kn+Ndc+NY 


(21) 


where  A^signai  1®  the  number  of  signal  electrons  in  Eq.  (18), 
ApN  is  the  number  of  photon  noise  electrons  in  Eq.  (19), 
Ndc  is  the  number  dark  current  noise  electrons  in  Eq.  (20), 
and  Nk  is  the  number  of  readout  noise  electrons  specihed 
by  a  manufacturer. 

The  SNR  can  be  dehned  again  expanding  all  the  equa¬ 
tions  for  noise  electrons: 


SNR= 


/|;2£(\)d\[A(D2/D,y](1.22\/4)2[To(\//2c)QE]7’,. 
{/;;2£(\)d\[A(D2/D2f,)](1.22\/4)2[ro(X//2c)QE]r,.  +  (RDcE,)+A2}i/2- 


(22) 


Using  Eq.  (22),  the  optimal  frame  rates  and  sensitivities  can 
be  obtained  to  achieve  the  required  SNR  for  a  given  target. 

2.3.5  Structural  misalignment  toierancing 

Using  the  mode  of  nonsequential  ray  tracing  of  ZEMAX,  a 
complete  ARGOS  optics  layout  is  constructed  based  on  the 
optical  speciheations  of  a  subaperture,  pyramidal  mirror, 
and  the  beam  combining  telescope,  as  shown  in  Pig.  11.  We 
intentionally  perturb  the  subtelescope  or  pyramidal  mirror 


to  determine  allowable  structural  misalignment,  and  we 
compensate  the  tilt  error  by  changing  the  tilt  angle  of  the 
fold  mirror  attached  to  the  PSM.  At  0.01 -deg  tilt  of  a  sub¬ 
aperture,  a  pure  PSM  motion  cannot  restore  the  SR  above 
0.8.  But  the  addition  of  PSM  piston  motion  can  restore  the 
SR  value  to  0.859.  We  could  achieve  a  SR  of  0.859  (which 
is  above  diffraction  limited)  over  a  0.01 -deg  tilt.  But  due  to 
a  magniheation  factor  10,  the  PSM  compensation  exceeded 
its  max  range  (0.6  mrad=  0.034  deg).  Since  we  mounted  a 
PSM  onto  a  precision  tip-tip  mount,  which  is  capable  of 
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Fig.  11  ZEMAX  3-D  nonsequential  ray  tracing. 


several  arcsecond  adjustment  (see  Table  2),  this  static  error 
does  not  limit  the  FSM  performance.  However,  it  is  much 
safer  to  have  a  FSM  within  a  range  of  eliminating  a  pos¬ 
sible  maximum  alignment  error,  0.005  deg  or  15  arcsec  for 
subtelescope  structural  misalignment  is  suggested. 

By  assuming  that  all  other  optical  components  are  per¬ 
fectly  aligned  and  the  FSM  can  compensate  for  all  the  re¬ 
sidual  tilt  errors,  the  tilt  errors  for  each  surface  of  the  py¬ 
ramidal  mirror  are  calculated.  When  the  tilt  error  of  the 
pyramidal  mirror  unit  equals  the  tilt  compensation  of  a 
FSM,  the  aberration  loss  due  to  the  tilt  is  completely  elimi¬ 
nated.  Therefore  there  is  no  theoretical  tilt  tolerance  for  the 
pyramidal  mirror  as  long  as  it  does  not  exceeds  the  maxi¬ 
mum  compensation  range  (0.01  deg). 

The  beam  combiner  was  tilted  along  the  x  and  y  axes 
while  leaving  other  optical  components  perfectly  aligned 
(Fig.  12).  This  beam  combiner  misalignment  is  not  correct¬ 
able  by  optical  actuators  like  FSMs.  However,  it  turns  out 
that  we  can  tolerate  up  to  0.2  deg  for  the  beam  combiner, 
which  is  less  stringent  than  other  misalignment  tolerances. 

3  Control  System 

3.1  Attitude  Control  System  Overview 

The  hnal  ARGOS  structure  and  the  major  ACS  components 
are  depicted  in  Fig.  1.  The  FOV  of  the  main  CCD  is  3 
aremin  and  thus  to  give  us  1/2  aremin  of  margin  on  either 
side,  the  ACS  subsystem  is  required  to  provide  a  pointing 
accuracy  of  ±  1  aremin.  The  period  of  operation  of  the  AR¬ 
GOS  system  without  human  intervention  must  be  60  min  or 
greater  (see  Table  1),  meaning  that  the  ACS  system  will 
have  to  either  not  saturate  its  actuators,  or  have  some  way 
of  desaturating  them  within  this  given  time  span.  The  sys- 


Fig.  12  PSF  plots  when  the  beam  combiner  has  tilt  errors.  From  left 
to  right,  X  tilt:  0.2,  Y  tilt:  0.4;  X:  0.3  Y:  0.3;  X:  0.25  Y:  0.25  (deg). 


Range:  50  deg  Resolution:  1 .5  deg 

ACS  Sensor 


7  X  5.25  deg 
FOV 


Viewfinder 


f=38.8  mm 


3.07  X  2.3  aremin 
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Fig.  13  Multistaged  ACS  sensors. 


tern  must  be  able  to  slew  at  a  rate  of  at  least  1.5  deg/s, 
placing  a  minimum  requirement  on  the  capabilities  of  the 
reaction  wheels  to  slew  the  spacecraft.  Due  to  the  nature  of 
the  air  bearing  system  chosen  to  simulate  the  space-based 
operation  of  ARGOS,  the  center  of  gravity  and  the  center  of 
rotation  of  the  body  will  not  necessarily  be  at  the  same 
position.  Three  perpendicular  precision  linear  slides  of  the 
active  balancing  system*^  (ABS)  correct  for  small  center  of 
mass  offsets  that  would  hurt  the  closed-loop  performance. 
The  ABS  can  also  be  used  to  remove  momentum  from  the 
wheels  by  intentionally  causing  an  offset  in  the  appropriate 
direction  (see  Fig.  1). 

The  sensor  suite  is  composed  of  three  integral  elements. 
First,  the  TCM-2-50  electronic  compass,  with  a  three-axis 
magnetometer,  a  two-axis  tilt  sensor,  and  facilities  to  pro¬ 
vide  temperature  information.  Inclinometers/electronic 
compasses  measure  the  relative  angles  between  the  inertial 
coordinate  frame  and  the  body  hxed  frame.  That  is,  they  are 
used  to  give  relative  elevation,  roll,  and  azimuth  informa¬ 
tion  between  these  two  coordinate  systems.  The 
inclinometer/electronic  compass  is  an  essential  component 
of  the  coarse  pointing  sensor  suite.  Second,  there  is  an  in¬ 
termediary  sensor,  which  takes  the  form  of  a  scope.  It  is  the 
intermediate  wide  FOV  CCD  that  provides  sufficient  over¬ 
lap  (greater  than  4)  between  the  TCM  (Fig.  13)  and  the 
main  imager.  Third,  there  is  a  three-axis  rate  gyroscope  for 
angular  rate  measurement. 

3.2  Control  Hardware  and  Avionics 

The  HEPC8  digital  signal  processing  module  houses  a  TI 
C6701  digital  signal  processor  (DSP),  two  12-bit  eight- 
channel  digital-to-analog  (D/A)  (three  signals  for  each 
FSM  and  three  RWA)  and  one  14-bit  eight-channel  analog 
to  digital  (A/D).  Each  reaction  wheel  speed  command  from 
the  DSP  is  sent  through  smoothing  hlters  with  a  cutoff  fre¬ 
quency  of  280  Hz.  The  internal  analog  motor  controller 
employs  a  proportional  integral  derivative  (PID)  logic  us¬ 
ing  the  motor  encoder  sensor  outputs.  A  set  of  three  piezo¬ 
electric  stacks  for  a  single  aperture  ESM  (a  total  of  nine 
piezolinear  stacks)  is  driven  by  the  D/A  signals,  which  are 
hltered  by  a  303-Hz  cutoff  frequency  low-pass  hlter  and 
then  amplihed  by  the  compact  voltage  ampliher  (PI 
E0669.OE)  located  in  the  right  side  of  ARGOS.  The  rate 
gyroscope  measures  the  rotational  rates  of  ARGOS  and 
feedbacks  the  voltages  to  the  A/D  with  a  resolution  of 
50  mV/(deg/s).  An  antialiasing  low-pass  hlter  with  a  48-Hz 
cutoff  frequency  was  implemented  to  reduce  the  high- 
frequency  noise  sent  to  the  A/D  channels.  The  electronic 
compass  is  directly  connected  to  the  PC  motherboard  using 
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WHEEL-3 


Fig.  14  System  identification  of  ARGOS-R/W  3  to  gyros  1 , 2,  and  3  outputs  (frequency  in  hertz). 


a  standard  RS232  serial  port.  The  PC -DSP  first  in  first  out 
(FIFO)  designated  for  direct  communication  between  the 
host  PC  memory  and  the  DSP  memory,  sends  the  image 
processing  data*^  such  as  the  centroids  of  the  target  from 
the  intermediate  CCD,  and  the  compass  readings  to  the 
DSP  controller  in  real  time.  To  meet  the  requirement  of 
simultaneous  automatic  control  of  various  subsystems,  the 
DSP  modules,  the  ABS  controller  and  the  wireless  card  for 
the  communication  with  a  remote  laptop  control  ground 
station,  all  reside  in  the  peripheral  component  interconnec¬ 
tion  (PCI)  slots  of  the  Athlon  1.4-GHz  PC  system  with 
1024  MB  RAM.  The  graphic  user  interface  (GUI)  PC  con¬ 
trol  application  sends  the  reference  slew  speed  and  pointing 
(azimuth  and  elevation)  for  a  specific  target  to  the  DSP 
while  performing  real-time  image  processing  from  both  the 
intermediate  ACS  CCD  and  the  main  imager  (spot-based 
wavefront  sensing^^). 

3.3  Pointing  and  Tracking  Controi  Design 

The  objective  of  system  identification  is  to  determine  the 
dynamic  properties  of  ARGOS  and  correctly  design  an  ap¬ 
propriate  control  system  to  capture  all  present  issues  that 
would  deteriorate  the  performance.  The  3X3  transfer  func¬ 
tion  matrix  between  the  three  D/A  signals  from  the  DSP 
and  the  three  gyro  outputs  to  the  A/D  port  for  the  reaction 
wheel  assembly  (RWA)  control  is  obtained.  In  this  case,  the 
system  identification  will  give  us  a  transfer  function  with 
the  ACS  filter  board  reference  command  as  the  input,  the 
smoothing  low-pass  filter,  motor  controller,  reaction 
wheels,  structure,  rate  gyros,  and  antialiasing  low-pass  filter 
as  the  plant,  and  the  filtered  gyro  signal  from  the  ACS  filter 
board  as  the  output. 

The  transfer  function  of  the  plant  is  then  found  by  di¬ 
viding  the  cross  spectral  density  of  the  input  and  output  by 


the  power  spectral  density  of  the  input.  Then,  the 
measurement-based  system  identification  algorithm  called 
integrated  frequency  domain  observability  range  space  ex¬ 
traction  and  Least  Square  parameter  estimation 
algorithm*®'*^  (IFORSELS),  is  performed  to  derive  a  state- 
space  representation  of  the  system  dynamics.  IFORSELS 
algorithm  integrates  the  frequency  domain  observability 
range  space  extraction  (EORSE)  identification  algorithm, 
the  balanced  realization  (BR)  model  reduction  algorithm, 
and  the  logarithmic  and  additive  least  squares  (LLS)  modal 
parameter  estimation  algorithms  for  low-order  highly  accu¬ 
rate  model  identification. 

The  EORSE  algorithm  is  a  slight  variation  of  the  sub¬ 
space  identification  algorithm  derived  by  De  Moor  et  al.*® 
and  its  objective  is  to  minimize  the  following  cost  to  obtain 
the  system  matrices.  A,  B,  C,  D, 

K 

7=2  \\G{io,)-[C{cxp{jw,At)I-A)]-^B  +  D)\\i  (23) 

k=l 

where  G{u)i^  is  the  frequency  response  samples  from  ex¬ 
periments. 

The  BR  model  reduction^'*  algorithm  transforms  a  state 
space  model  to  a  balanced  coordinate  and  reduces  the 
model  by  deleting  the  states  associated  with  the  smallest 
balanced  singular  values.  The  LLS  estimation  algorithms 
improve  the  fitting  of  reduced  models  to  experimental  data 
by  updating  state  space  parameters  in  modal  coordinates. 

Erom  the  experimental  transfer  function  measurements, 
a  0.01-  to  40-Hz  chirp  signal  generated  the  sampled  transfer 
functions  in  Fig.  14  (solid  line).  Using  the  FORSE  algo¬ 
rithm,  an  80th-order  state-space  model  was  selected  to  ap¬ 
proximate  the  experimental  data  and  then  BR  and  LLS  al- 
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gorithms  were  used  to  get  a  60th-order  IFORSELS  model 
(shown  by  dotted  lines  in  Fig.  14).  The  low-frequency 
peaks  below  0.1  Hz  is  the  pendulum  mode  of  ARGOS 
while  the  first  lightly  damped  structural  mode  occurs  at  8  to 
9  Hz.  Based  on  the  linear  model  obtained,  the  pointing  and 
tracking  controller  was  successfully  designed  and  imple¬ 
mented  performing  a  large  angle  slew  maneuver.^*  The 
controller  calculates  the  error  quaternions^^  from  the  rate 
gyroscope,  the  electronic  compass/inclinometer  and  the 
viewhnder  centroiding  CCD  mixing  the  mechanical  sensors 
with  the  optical  sensor  (the  viewhnder  CCD  in  Figs.  1  and 
13). 

4  Conclusion 

The  ARGOS  is  successfully  integrated  into  the  full  struc¬ 
ture  ready  to  operate.  The  ARGOS  testbed  is  the  hrst 
sparse-aperture  array  simulating  a  space-borne  observatory 
in  a  1-g  environment  representing  real-world  problems 
such  as  the  vibrational  coupling  between  a  spacecraft  struc¬ 
ture  and  the  wavefront  errors  propagating  through  the 
whole  system.  In  this  paper,  the  physical  nature  of  inter¬ 
ferometric  arrays  was  studied  and  the  design  trades  to  meet 
the  performance  requirements  were  explored.  A  preliminary 
assessment  showed  that  the  beam  combining  problem  is  the 
most  challenging  aspect  of  sparse  optical  arrays.  The  need 
for  optical  control  is  paramount  due  to  tight  beam  combin¬ 
ing  error  tolerances.  The  attitude  controller  design  in  this 
paper  is  focused  on  the  attitude  control  system  incorporat¬ 
ing  the  coarse  optical  sensor  and  general  ACS  sensor  suite. 
The  interferometric  beam-combining  controller  utilizing 
model-based  control  and  neural  network  is  under  develop¬ 
ment  as  well  as  software-based  wavefront  sensing  tech¬ 
niques.  The  wavefront  sensing/control  requirements  appear 
to  be  a  major  technology  and  cost  driver.  We  anticipate  that 
ARGOS  will  provide  an  experimental  platform  for  the  de¬ 
velopment  of  many  phasing  and  vibration  control  tech¬ 
niques  for  future  interferometers.  Follow-on  work  will  fo¬ 
cus  on  actual  versus  expected  imaging  performance  for  a 
range  of  targets. 
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